Preprint typeset in JHEP stylo - HYPER. VERSION 

Superradiant instabilities of rotating black branes and strings 



in ■ 

o : 

O ! Vitor Cardoso^ >* and Shijun Yoshida ft 

CM. 

'McDonnell Center for the Space Sciences, Department of Physics, Washington University, St. Louis, 
Missouri 63130, USA 

* Centro de Fisica Computacional, Universidade de Coimbra, P-300^-516 Coimbra, Portugal 
''Science and Engineering, Waseda University, Okubo, Shinjuku, Tokyo 169-8555, Japan 



X 



vcardosoOwugrav.wustl . edu, 3hijun@waseda. jp 

(N 
(N 
O 

Abstract: Black branes and strings are generally unstable against a certain sector of gravitational pertur- 
bations. This is known as the Gregory-Laflamme instability. It has been recently argued || that there 
42 1 exists another general instability affecting many rotating extended black objects. This instability is in a 
sense universal, in that it is triggered by any massless field, and not just gravitational perturbations. Here 
we investigate this novel mechanism in detail. For this instability to work, two ingredients are necessary: 
42 ■ (i) an ergo-region, which gives rise to superradiant amplification of waves, and (ii) "bound" states in the 
\ effective potential governing the evolution of the particular mode under study. We show that the black 
brane Kerr4 x R p is unstable against this mechanism, and we present numerical results for instability 
timescales for this case. On the other hand, and quite surprisingly, black branes of the form Kerr^ x R p 
are all stable against this mechanism for d > 4. This is quite an unexpected result, and it stems from 
the fact that there are no stable circular orbits in higher dimensional black hole spacetimes, or in a wave 
picture, that there are no bound states in the effective potential. We also show that it is quite easy to 
simulate this instability in the laboratory with acoustic black branes. 
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1. Introduction 

We are all painfully aware that exact, closed form solutions to physical problems are an exception. The 
rule is that one must resort to all kinds of tricks and assumptions in order to get a grip on the equations: 
this is why, despite their little insight into the physics, numerical solutions are proliferating as problems 
become more and more complex. Exact solutions are therefore most welcome. Exact solutions to Einstein 
equations are extremely useful, especially if they describe simple yet physically attainable systems. Indeed, 
take for example the famous Schwarzschild metric: with this exact solution at hand, describing the geometry 
outside a spherically symmetric distribution of matter, one was able to compute the deflection of light as 
it passes near the Sun (and to match the theoretical prediction against the observational data), thereby 
giving strong support to Einstein's theory. We now know that the outside geometry of many astrophysical 
objects is well described by the Schwarzschild metric, and we can start studying them by investigating the 
properties of this metric. 

One of the most important things that one should study first is the classical stability of a given 
solution: if a solution is unstable, then it most certainly will not be found in nature (unless the instability 
is secular) and the solution loses most of its power. What does one mean by stability? In this classical 
context, stability means that a given initially bounded perturbation of the spacetime remains bounded for 
all times. For example, the Schwarzschild spacetime is stable against all kinds of perturbations, massive or 
massless §. Thus, if one considers small deviations from the Schwarzschild geometry (for instance throwing 
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a small stone into the black hole, this disturbs the geometry), stability means that these small perturbations 
eventually fade away. They will never disrupt the spacetime. The term "perturbations" means either small 
(metric) deviations in the geometry or small-amplitude fields in the geometry, for example scalar fields or 
electromagnetic fields (as such, deviations in the geometry can be looked at as perturbations induced by 
gravitons). So, the Schwarzschild geometry is indeed appropriate to study astrophysical objects. 

Here we shall investigate an instability which is rotation-triggered, and which was hinted at for the 
first time by Marolf and Palmer The mechanism has recently been explained in ||], where it was 
shown that it could be triggered in many systems. However, a detailed description is still lacking, and 
this is the main purpose of the present work. The physical nature of this instability can be understood 
with the following facts: It is known that the Kerr geometry displays superradiance This means that 
in a scattering experiment of a wave with frequency to < mVL (m is again the azimuthal wave quantum 
number) the scattered wave will have a larger amplitude than the incident wave, the excess energy being 
withdrawn from the object's rotational energy. Now suppose that one encloses the rotating black hole 
inside a spherical mirror. Any initial perturbation will get successively amplified near the black hole event 
horizon and reflected back at the mirror, thus creating an instability. This is the black hole bomb, as 
devised in || and recently improved in ||. This instability is caused by the mirror, which is an artificial 
wall, but one can devise natural mirrors if one considers massive fields J?], ||. Imagine a wavepacket of the 
massive field in a distant circular orbit. The gravitational force binds the field and keeps it from escaping 
or radiating away to infinity. But at the event horizon some of the field goes down the black hole, and if the 
frequency of the field is in the superradiant region then the field is amplified. Hence the field is amplified 
at the event horizon while being bound away from infinity. Yet another way to understand this, is to think 
in terms of wave propagation in an effective potential. If the effective potential has a well, then waves get 
"trapped" in the well and amplified by superradiance, thus triggering an instability. In the case of massive 
fields on a (four-dimensional) Kerr background, the effective potential indeed has a well. Consequently, 
the massive field grows exponentially and is unstable. It is the presence of a bound state that simulates 
the mirror, and without a bound state we should never get an instability. 

Now, consider for example the following rotating black-brane geometry 

ds 2 = ds^err + dx % dXi , (1-1) 

where Kerr stands for the usual Kerr geometry. It's a known property, and we will show it explicitly, that 
the propagation of a massless field (scalar, electromagnetic, or gravitational) in this black-brane geometry, 
is equivalent to the propagation of a massive field in the vicinity of the Kerr black hole. Thus, the particular 



black brane (IT) is unstable 



The instability argument applies to many rotating extended objects, and here we shall study some of 



them, making an extensive analysis of the d— dimensional case (IT). We show that for d = 4 they are 



unstable, and we present detailed results on the stability. For d > 4 the black branes described by 

ds 2 = ds 2 ^ erid + dx l dxi , (1-2) 

are stable, where Kerrd is the Myers-Perry || rotating black hole. This is due to the non-existence of 
stable bound orbits for massive particles or, in terms of wave propagation, there is no well in the effective 
potential for these systems. To conclude, we show that one can mimic this instability in the laboratory 
using analogue acoustic black branes. 
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2. Rotating Kerr-like black branes 



In this section we study in detail the black branes of the form (|T]^) . In higher dimensions there are several 
choices for rotation axes of the Myers- Perry solution (labeled Kerrj in (|1.2j)) and there is a multitude of 
angular momentum parameters, each referring to a particular rotation axis We shall concentrate on 
the simplest case, for which there is only one angular momentum parameter, denoted by a. 

2.1 Formalism 

In Boyer-Lindquist-type coordinates the black branes we study in this section are described by 

,2 



A -a 2 sin 2 9 2a(r 2 + a 2 - A) sin 2 
ds z = dt 2 i dtdtp + 



,2 + a 2 )2 _ Aa 2 sin 2A v , fljn2 



sin^ 6dtp* + -^dr z + Ed0 z + r z cos z Od^i + dx L dx u (2.1) 



where 



E = r 2 + a 2 cos 2 (9, A = r 2 + a 2 - Mr 1 "", (2.2) 

and dfi 2 denotes the standard metric of the unit n-sphere (n = d — 4), the x l are the coordinates of 
the compact dimensions, and % runs from 1 to p. This metric describes a rotating black brane in an 
asymptotically flat, vacuum space-time with mass and angular momentum proportional to M and Ma, 
respectively. Hereafter, M, a > are assumed. 

The event horizon, homeomorphic to S 2+n , is located at r = r + , such that A| r=r+ = 0. For n = 0, an 
event horizon exists only for a < M/2. When n = 1, an event horizon exists only when a < \[M , and the 
event horizon shrinks to zero-area in the extreme limit a — * \f~M . On the other hand, when n > 2 , A = 
has exactly one positive root for arbitrary a > 0. This means there is no bound on a, and thus there are 
no extreme Kerr black branes in higher dimensions. 

Consider now the evolution of a massless scalar field ^ in the background described by ([2.1|), The 
evolution is governed by the curved space Klein-Gordon equation 

° 1 ^r4rA=^ (2-3) 



dx^ \ dx u 

where g is the determinant of the metric. We can simplify considerably equation ( |2.3D if we separate the 
angular variables from the radial and time variables, as is done in four dimensions This separation 
was accomplished, for higher dimensions, in [11] for five dimensional Kerr holes and also in [|l^] for a 
general 4 + n-dimensional Kerr hole. Since we are considering only one angular momentum parameter, the 
separation is somewhat simplified, and we can follow [13|. In the end our results agree with the results in 
[11, [l^] with only one angular momentum parameter in their equations. 

We consider the ansatz <f> = e~ lu)t+imv+1 ^ iXt R(r) S (9)Y (tt) , and substitute this form in ( |2~3|) , where 
Y(p,) are hyperspherical harmonics on the n-sphere, with eigenvalues given by —j(j+n—l) (j = 0, 1, 2, • • • ). 
Then we obtain the separated equations 



1 ( d . „ „ n dS 
sm 



sin cos n \d9 dB 
and 



\ + [a 2 (cu 2 - fi 2 ) cos 2 - m 2 esc 2 - j(j + n - 1) sec 2 + A] S = 0, (2.4) 



t d f nA d_R\ + f nr 2 + a 2 )-ma] _ j(j + n -l)a 2 _ ^ _ 
dr V dr J A r z 
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where A is a constant of separation, A := A — 2muja + to 2 a 2 , and /x 2 = Y^f^i- Interestingly, note the 
important point that Equations ( |2.4j )-( |iE5D are just those that describe the evolution of a massive scalar 
field, with mass /i, in a d— dimensional Kerr geometry. 

The equations (|2.4| ) and (|2.5| ) must be supplemented by appropriate boundary conditions, which are 
given by 

(r — r_|_) _lcr as r — > r + 

as r — > oo , 



with 



[(r 2 + a 2 )w - ma] r + 

(T = 7Z — T^r^^—^JT ■ ( 2 - 7 ) 



n 



l)(r 2 +a 2 )+2r 2 



In other words, the waves must be purely ingoing at the horizon and purely outgoing at the infinity. For 
assigned values of the rotational parameter a and of the angular indices I ,j ,m there is a discrete (and 
infinite) set of frequencies called quasinormal frequencies, QN frequencies or cjqjv, satisfying the wave 
equation ( f2~5| ) with the boundary conditions just specified by Eq. Q2.6|) . Since the time-dependence of 
the mode is e~ lut , unstable modes will have frequencies with a positive imaginary part, and thus grow 
exponentially with time. 

2.2 The physical nature of the instability 

The equations presented in the previous section are amenable to numerical calculations (infra), but don't 
offer much physical insight into the physics. To gain some more intuition, and understand why we expected 
an instability to be triggered, we find it convenient to make instead the following change of variables, 
following Furuhashi and Nambu ||: 

R = r -n/2-ly^ ^ 



in which case the wave equation (|2.q ) is transformed into 

d 2 ^ 



dzl 



V^ = 0. (2.9) 



Here, we have defined the modified tortoise coordinate as J~ = -5. (Note that the usual definition of 
the tortoise coordinate r* is dr/dr* = A/(r 2 + a 2 )). 
The effective potential V is equal to 

I Tt 71/ (X^ Tip 1 

-V = -tt(-X - a 2 [i 2 2amuj + 2a 2 w 2 ) + — r(a 2 a; 2 - 2amuo jn + m 2 - A - j 2 + j + 2) + 

r z 2 4 r 4 2 

a 4 , 9 n n 2 u 2 M M n N M 2 . n 2 , 

-(2 + j-j 2 + --nj--) + ^ + ^{\-l--) + — i (l + n+-) + 

Ma 2 , n 3n, o , 

^ T 50- 2 -j-3 + jn + T ) + W 2 - / , 2 . (2.10) 

Upon close inspection, this potential has two key ingredients which can trigger the instability: super- 
radiant amplification of scattered waves and the existence of bound states. We shall briefly explore and 
describe these features in what follows. 



-A- 



2.2.1 Superradiant scattering 

In a scattering experiment (with u) > fj, so that plane waves at infinity are possible), equation ( |2.9| ) has the 
following asymptotic behavior: 

{T(r — r+)~ lcr as r — > r+ , , 

• AT^2 ■ AI2 — 2 2.11 

where a was defined in the previous section, equation ([2.7|). 

These boundary conditions correspond to an incident wave of unit amplitude from +00 giving rise to 
a reflected wave of amplitude R going back to +00 and a transmitted wave of amplitude T at — 00 (as 
before, the boundary conditions we impose do not allow for waves emerging from the horizon). Since the 
potential is real, the complex conjugate of the solution ^1 satisfying the boundary conditions (|2.11[) will 
satisfy the complex-conjugate boundary conditions: 

\T*(r — r + ) 1 ' 7 as r — > r + , 

*2~< ■ HT~ — 2 . 2 (2-12) 

i?*e _1 v w v r + e l v w v- 



r as r —>■ 00 . 



Now, these two solutions are linearly independent, and standard theory of ODE tells us that their Wron- 
skian is a constant (independent of r). If we evaluate the Wronskian near the horizon, we get 

W = -2iaA\T\ 2 (2.13) 

where the constant A is defined by 

A ,r^r+, (2.14) 



r 2 (r — r+) 

and is positive definite. Evaluating the Wronskian at infinity we get 

W = 2iV^ 2 - V 2 (\R\ 2 ~ 1) • (2.15) 

Equating the two we get 



|i?| 2 = 1 \T\ 2 . (2.16) 

Vw 2 -// 2 



Now, in general \R\ is less than unity, as is to be expected. However, for 

co-^^<0, (2.17) 

we have that a is negative, and therefore in this regime (which will be referred to as the superradiant 
regime) |i?| 2 > 1. This means that one gets back more than one threw in. Of course the excess energy 
comes from the hole's rotational energy, which therefore decreases. As a final remark, we notice that we 
have been assuming without loss of generality uj > 0, and thus superradiance only exists for m > 0, which 
are modes co-rotating with the black hole. 
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Figure 1: A typical form for the effective potential, here shown for I = m = 1 modes. We have set M = 1, so 
the rotation parameter a varies between (Schwarzschild limit) and 1/2 (extremal limit). Here we plot the effective 
potential for the near extreme situation, a ~ 0.5 and for /i = 0.7 and ui = 0.6878. 



2.2.2 The potential well 



The effective potential in four dimensions, d = 4, as given by ( 2.10| ) is plotted in figure [l], for a ~ 0.5 



(in units of M), fj, = 0.7 and uj = 0.6878. As can be seen from Figure [l] the potential has, in the four- 
dimensional situation, two extremum between the event horizon and spatial infinity. The local minimum 
creates a "well" -like structure, which will be so important to trigger the instability. The potential is 
asymptotic to (fi 2 - co 2 ) at spatial infinity. That a well must necessarily arise in four dimensions can be 



seen from the asymptotic nature of the potential. In fact, for n = and for large r, the potential in ( 2.10 ) 
behaves as 

V ~/i 2 -uj 2 -ii 2 Mr~ x , r -» oo , n = 0. (2.18) 

Thus, the derivative is V ~ ^-j^, which is positive. Thus a well necessarily arises. 

The instability can now be easily explained: the waves get "trapped" inside the potential well and 
amplified by superradiance. Thus, an exponential growth of any initial perturbation is inevitable, and 
therefore an instability arises. This simple picture is accurate enough that it will allow us to predict if and 
when an instability is triggered. 



2.3 The numerical search technique 

The problem we have to solve numerically is a simple boundary value problem of a second-order ordinary 
differential equation. However, this is not easy to solve with the usual techniques of direct integration, 
because of difficulties related to the absence of incoming radiation at the boundaries flj]. A standard 
technique to overcome the difficulties is the so-called Leaver's method H], which we employ in this paper. 
In fact, this seems to be the first attempt at trying to find instabilities using this method, and it works 
well, as we shall see. In what follows, we shall describe the four dimensional (n = 0) case thoroughly, for 
concreteness. The method is easily adapted to any number of dimensions n, but the equations get much 
lengthier, so we refrain from presenting them here. We pick units such that M = 1. First, let us consider 
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the radial equation. Following Leaver, the radial function R can be expanded around the horizon as 



oo , _ \ k 

R=(r- r + )- la (r - r _)-i+^+^ 2 (2^i)- 1 +^i e ^ 1 r £ ftfe I r_r+ \ ^ 

k=o \ r ~ r -J 



where a and u>i are defined by 



(r 2 , + a 2 )co — ma /—^ = 

° = + 2 '2 r +i "1 = V^-fi 2 . (2.20) 



Here, 60 is taken to be 60 = 1> and r + and r_ are, respectively, the coordinate radius of the event and the 
Cauchy horizon of the Kerr black hole, given by r± = (1 ± b)/2, where b = (1 - 4a 2 ) 1 / 2 . Note that the 
branch of \fz has been chosen such that —ir/2 < arg^/z < tt/2. The expansion coefficients bk in equation 
( |2.19| ) are determined via the five-term recurrence relation, given by 



where 



a h + p b = , a\b 2 + fi\b\ + jib = , a 2 ^3 + /?2&2 + 72&1 + S 2 bo = , 

OLkbk+i + Pkbk + lkbk~i + hbk-2 + = , k = 3, 4, 5, ■ ■ ■ , (2.21) 



a fc = (1 + fc)(l + fc + /i ) , A = -4A; 2 + + oi , 7 fc = 6A; 2 + /i 2 /c + o 2 , 
5fe = -4/c 2 + /i 3 £; + 03 , efc = k 2 + /i 4 A; + 04 . 



(2.22) 



Here, the coefficients appearing in equation ( |2.22| ) are given by 
h = -ib~ l {-2am + oj(1 + b)} , 

h x =-2- Siamb- 1 + i{4u;6 _1 (l + b) + ^Ajjf 1 + 2u>i(l + 6)} , 

/i 2 = -i(6wi) _1 [6wi(a; - 2am) + 46 2 w 2 + Sfoj^ 2 + 2wj(w + a>i - »)}] , (2.23) 
/i 3 = 10 - %iamb~ l + i{4w(l + o" 1 ) + 3//^ 1 + 2u>i(3 + b)} , 
«4 = — i{buji)~ l [u)i(u) — 2am) + b{[i 2 + uj\(oj + 2oj\ — 4i)}] , 
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and 

01 = (4bujl)- 1 [b 3 ujf + 2(2am-uj){-2uj 2 {i + uj + uj 1 ) + fi 2 (2i + 2uj + ui)} 

-2b 2 uj\{ii 2 - 2(u? + + w))} 

+b{fi 4 + /u 2 (4 + 4j4; m — Aiuj — 8uj 2 — 2iu)\ + 8amuJi — Qujui) 
+ lo 2 {-A - 4A lm + 4ilo + 7lo 2 + Aiuox - 8amwi + 8wo>i)}] , 

2 = (4bu%)~ 1 [-2l?u( - 46 2 {/u 4 - fi 2 {3uj 2 + 2iwi + 2wu>i) + 2(u 4 + + u 3 oji - 22w 3 )} 

-6(2am - w){/i 2 (2z + 2a; + Wi) - 2(w 3 - 2iuj\ + w 2 (z + 

+b{— 3/i 4 — 2// 2 (6 + 4A/ m — 6iw — 10u> 2 — lOiwi + 8amwi — 7uui) 

+ 2(-9w 4 + \2iujuj\ + 12iw 3 - 2w 3 (3i + 5wi) + w 2 (6 + 4A im - 8wi + 8amwi))}] , (2.24) 

3 = (46w 2 )~ 1 [6 3 w 4 + b{3fj, 4 + 15w 4 - 32iuw 2 - 16w 2 (l + &i) + 4w 3 (3i + 4u>i) 

— 4w 2 (3 + Ai m — hiuj\ + 2amwi) + 2 / u 2 (6 + 2^ m — 6iu — 8uj 2 — YJiui + Aamuji — buui)} 
+26 2 {/i 4 - ^ 2 (3w 2 + 2iux + 2uux) + 2(w 4 + iu?0Ji + u^ujx - Aicof )} 
+2(2am - w){3^ 2 (2i + 2w + u/i) - 2(3w 3 - 8iw 2 + 3uj 2 (i + wi))}] , 

4 = (46w 2 ) _1 [-2(2am-a;){/i 2 (2i + 2a; + wi) - 2(w 3 - 3ia; 2 + w 2 (i + ui))} 

+6{-/i 4 - 4(w 4 + oj 2 {-1 + 2iui) - 3iuuj\ - 3w 2 (l + 2iu x ) +co 3 {i + ui)) 
+2fi 2 (-2 + 2u? + 8icji + w(2i + wi))}] , 

Using Gaussian elimination twice, one can reduce the five-term recurrence relations ( |2.21| ) to three-term 
recurrence relations, which can be written as 

«o ft i + $) & o = , a' k b k+ i + (3' k b k + 7^-1 , k = 1, 2, • • • , (2.25) 

Since the Gaussian elimination is straightforward, we do not show an explicit procedure for the Gaussian 
elimination (see, e.g., [15]). An eigenfunction satisfying the quasinormal mode boundary conditions behaves 
at the event horizon and infinity as 



(r — r+) la as r — ► r 



* ~ I ^-l-H/5(2a)i)- 1 +K.)i e iwir as r ^ o^. ( 2 - 26 ^ 

Therefore, one can see that the expanded wave function (|2.19| ) satisfies the quasinormal mode boundary 
conditions if the expansion in (2.19) converges at spatial infinity. This convergence condition for the 
expansion ( |2.19[), namely the quasinormal mode conditions, can be written in terms of the continued 
fraction as p6| , [14}| 

, q' o7 ; ol x i 2 a' 2 ^' 3 _ , »o7i _ n (0 07 , 



where the first equality is a notational definition commonly used in the literature for infinite continued 
fractions. Here, we shall adopt this convention. 
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Next, we turn ourself to the angular equation. In order to determine the separation constant Ai m , a 
similar technique to that used for the radial equation can be applied. Imposing regularity on the symmetry 
axis, we can expand the angular function as 



S = e a ^ u (l-u 2 )^ i J2d k (l + u) k , (2.28) 



fc=0 



where u = cos 9, and do = 1. The expansion coefficients d k in equation ( 2.28 ) are determined via the 
three-term recurrence relation, given by 



where 



a di + o d o = , a k d k+1 + f3 k d k + %d k -i = , k = 1, 2, 3, • • • , (2.29) 
a k = -2(1 + k)(k + \m\ + 1), 

f3 k = k{k - 1) + 2k(\m\ + 1 - 2awi) - (2atJi - |m|)(|m| + 1) - {a 2 oj\ + A im ) , (2.30) 
7^ = 2auj\(k + \m\) . 

Note that in the four dimensional case, we have to set j = in equations (2.5) and (2.6). If the expansion 
in equation ( 2.28j ) converges for — 1 < u < 1, the angular function satisfies the regularity conditions at 



u 



±1. This condition is equivalent to the continued fraction equation, given by [JO], Q 



Pi- ft- Ps- 

Now that we have two continued fraction equations of the frequency u and the separation constant Ai m 
(ET27D and ( gjD , we can obtain the frequency u and the separation constant A[ m solving equations (|2.27|) 



and ( |2.3lD simultaneously. These coupled algebraic equations ( |2.27| ) and ( p.31| ) can be solved numerically 
(see, e.g., 0, H). 

2.4 The instability for the d = 4 case 

Some numerical results for d = 4 were given by Furuhashi and Nambu || and also by Strafuss and Khanna 
[||], in the context of massive field instability of the Kerr metric which as we saw, translate immediately 
to our situation. Nevertheless, we have performed an independent extensive numerical search, which 
allow us to have a more complete picture and understanding of the situation. Our numerical results were 
obtained using the technique just described supra, and are summarized in Figures |2[]3|. There is more 
than one frequency with a positive imaginary part, corresponding to unstable modes. We only show the 
most relevant one which is the one having a larger imaginary part, and therefore corresponds to the most 
unstable mode. 

In Figure [2] we plot the imaginary part of the (unstable) mode as a function of a for several values of 
H and for I = m = 1. The instability timescale is given by t ~ im [ cJ \ • ^ e can see that: 

(i) The instability gets stronger (the typical timescale decreases) as the rotation a increases. This is 
expected, since we know that the superradiant amplification gets stronger as a increases P, 19]. 



(ii) The instability also gets stronger as \x increases. This has a very simple explanation: the depth of 
the potential well in Figure |l| is larger for larger values of /i, and therefore the waves get more efficiently 
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Figure 2: Details for the instability in the four dimensional case. We have set M = 1, so the rotation parameter 
a varies between (Schwarzschild limit) and 1/2 (extremal limit). In the left panel we plot the imaginary part of 
the unstable modes as a function of a for nine values of [a, and for I = m = 1. The same quantities are plotted in 
the right panel, but now we focus on the maximum instability region, for values of /i near unity. The instability is 
stronger for larger values of a and for fi ~ 1. In fact, our results suggest that the imaginary part of the unstable 
modes attains its maximum value, 10 -6 x Re[w] for // = 0.9 and a — 0.497. Notice finally that as one decreases a 
(going right on the x axis in the left panel), there is a critical a below which there is no instability. Although there 
are an infinity of unstable modes, we only show here the most unstable one. 



trapped. Our results suggest that the imaginary part of the unstable modes attains its maximum value, 
1(T 6 x Re[w] for I = m = 1, n = 0.9 and a = 0.497. 

(iii) For a fixed \x and as one decreases a the unstable modes disappear below a certain critical value 
of a. This is because the superradiance condition Re[u;] — a 2™°s < is not satisfied for very low rotation 

parameter a. This is best seen in Figure ||. 

In Figure |3| we plot the real part of the (unstable) mode as a function of a for several values of \x and 
for I = m = 1. Two important features are the following: 

(i) As shown in both plots of Figure ||, the superradiant factor Re[w] — Jt^i is negative, and one is 

therefore in the superradiant regime. Thus, this is indeed a superradiant instability. 

(ii) The instability occurs only for u < fi. Otherwise, the waves would not be trapped in the potential 
well, but they would rather escape to infinity, and the perturbation would be damped. 

Detweiler |H| finds (in four dimensions), in the limit of small fj,, that the characteristic frequencies of 
the unstable modes are, in our units, 

a (uMf 

a; ~ « + iu— — — — . (2.32) 
p P M 3072 v ; 

Our numerical results fit this prediction very well. For near-extremal black holes, and /x ~ 0.1 Detweiler 's 

formula ( 2.32j ) predicts uj ~ 0.1 + 10~ 13 i, which checks very well with our numerics (see Figure | and |). 



Although of no direct interest for this work, we found also stable modes. Our results for stable modes 
will be published elsewhere ||23||, and compared with the existing analytical ones 



2.5 No instability for higher dimensional rotating Kerr-like black branes 

We have also performed a numerical search for the instability in the d = 5 , 6 , 7 case. We found no trace 
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Figure 3: More details of the instability in the four dimensional case. Here we plot the superradiant factor 
Re [a*] — a2 ™i as a function of a for I = m = 1 and for several values of /i. Again, we have set M = 1 and so 

the superradiant factor is, in four dimensions, Re [a;] — the instability disappears for a below a certain critical 
point, which is given by setting the superradiant factor to zero. On the right panel we make the same plot but 
on a logarithmic scale, for better visualization. Notice that the superradiant factor is negative, as it should for the 
instability to be triggered. 

of unstable modes, and we justify this in a rigorous analysis presented in Appendix |A|. We now turn to 
explain why higher dimensional Kerr-like branes should be stable against this particular mechanism. 

We saw that the instability arises because of superradiantly amplified trapped modes, in the potential 
well. But does a well exist for general d? It doesn't, and to understand this, we have to look at the 
asymptotic behavior of the effective potential ( 2,1C| ). If the derivative of this potential is positive near 



spatial infinity, we are guaranteed to have a well, and therefore an instability. If the derivative is negative, 
the modes should all be stable. Near infinity, the dominant terms in the effective potential fl2.10| ) are 

V ~ -\{A lm + aV 2 + 5 + ^- aV) + (n + l)fi 2 Mr- 2 - n , (2.33) 

where we have already substituted for the separation constant A = A[ m — 2amu> + a?uj 2 . It is immediately 
apparent that the four dimensional case is a special one: if n = 0, the dominant term in the derivative is 
(n + l)/i 2 Mr" 2 "", which is positive and we therefore are guaranteed to have a bound state. Thus, this 
case should be unstable, and it is, as we just described in the previous subsection. 

When n > the other terms dominate. In fact, for n > 1 they are positive. For n > 1 the dominant 
terms are 

V ~ -- 3 (A lm + aV + ^ + % - aV) . (2.34) 
r 6 2 4 

Since u < fj,, this is negative (the separation constant A[ m can be shown to be positive). Thus, for 
d > 5 there is no potential well, no bound states and therefore no instability, even though there is still 
superradiance. The situation for n = 1 is not as clear, because there is the extra term (n + l)/x 2 Mr -2 ~ n . 
In principle it should be possible to have, for certain very specific parameters, a potential well. But to do 
that, one would have to require that /i be very large, and this makes it very hard to study the problem 
numerically (the imaginary part is expected to be extremely small in this regime, as shown by Zouros and 
Eardley S, and this prevents any numerical treatment). 
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There is another, perhaps more physical, explanation of the fact that higher dimensional rotating 
black branes should be stable. We know that the a massless field propagating in the background of these 
rotating black branes is equivalent to a massive field propagating in the background of a Myers-Perry black 
hole. Imagine now a wavepacket of the massive field in a distant circular orbit. In the four-dimensional 
situation, the gravitational force binds the field and keeps it from escaping or radiating away to infinity. 
But at the event horizon some of the field goes down the black hole, and if the frequency of the field is in 
the superradiant region then the field is amplified. Hence the field is amplified at the event horizon while 
being bound away from infinity. We should thus have an instability, as we do. However, for d > 4, there 
are no stable orbits , and thus the field escapes to infinity, no instability is triggered. This complements 
the previous view point. Escaping to infinity is equivalent to not having a potential well to bound the field. 



3. Acoustic black branes 

In 1981 Unruh [25] introduced the notion of "dumb holes", which is an example of analogue black holes. 
While not carrying information about Einstein's equations, the analogue black holes devised by Unruh do 
have a very important feature that defines black holes: the existence of an event horizon. The basic idea 
behind these analogue acoustic black holes is very simple: consider a fluid moving with a space-dependent 
velocity, for example water flowing throw a variable-section tube. Suppose the water flows in the direction 
where the tube gets narrower. Then the fluid velocity increases downstream, and there will be a point 
where the fluid velocity exceeds the local sound velocity, in a certain frame. At this point, in that frame, 
we get the equivalent of an apparent horizon for sound waves. In fact, no (sonic) information generated 
downstream of this point can ever reach upstream (for the velocity of any perturbation is always directed 
downstream, as a simple velocity addition shows). This is the acoustic analogue of a black hole, or a dumb 
hole. We refer the reader to pq] for a review of these objects. 



3.1 The wave equation for a general space-dependent density 

We can also model acoustic black branes easily, and we can thus look for the instability dealt with here in 
the laboratory. For simplicity, let us consider the following black brane. First build a 2 + 1 dimensional 
acoustic black hole modelling a draining bathtub. Consider a fluid having (background) density p(x). 
Assume the fluid to be locally irrotational (vorticity free), barotropic and inviscid. From the equation of 
continuity, the radial component of the fluid velocity satisfies pv r ~ 1/r. Irrotationality implies that the 
tangential component of the velocity satisfies v e ~ 1/r. We shall keep both p and the sound velocity c as 
position-dependent quantities, thus generalizing the treatment in [26, |27{|. We then have 



= = (3.1) 

pr r 

where A , B are constants. The acoustic metric describing the propagation of sound waves in this "draining 
bathtub" fluid flow is IMj: 



/ A 2 1 o 2 + B 2 \ 2A 

ds 2 = - c 2 dt 2 + — drdt - 2Bd<j)dt + dr 2 + r 2 d<\> 2 . (3.2) 

\ r z J pr 

It is however better to work with a more transparent metric. Some physical properties of this draining 
bathtub metric are more apparent if we cast the metric in a Kerr-like form performing the following 
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coordinate transformation (see [27, pgj]): 



dt —>■ dt = dt =-= — f n / n dr (3.3) 

r 2 c 2 _ A 2 /p 2 v ' 

- ^ dr, (3.4) 



r(r 2 c 2 — A 2 / p 2 ) 
Then the effective metric takes the form 

d#* = -(l- A2,P l + B2 \ c 2 dt 2 +(l- 1 dr 2 - 2Bd4>dt + r 2 d4> 2 . (3.5) 

\ c z r z J \ p z c z r z J 



As explained in [27], this metric and the Kerr metric differ in an important aspect, in that whereas the 
rotation for the Kerr black hole is bounded from above, here it is not, at least in principle. Thus, B could 
be as large as desired. 



Now, we are free to add an extra dimension z to this effective geometry, as explained by Visser [26], 
resulting in the following black brane geometry 

ds 2 = -U- A2/ j£ B2 ) <? dP + ^272) 1 dr2 - 2Bd4>dt + r 2 d4> 2 + dz 2 . (3.6) 



The propagation of a sound wave in a barotropic inviscid fluid with irrotational flow, which is assumed to 
be the case, is described by the Klein-Gordon equation V^V^ = for a massless field f in a Lorentzian 
acoustic geometry. Separating variables by the substitution 

r , $) = > /r*(r)e < ^ +m *- w ^ , (3.7) 
implies that 'l'(r) obeys the wave equation 

d2 * + ((u-^-v)* = 0. (3.8) 



dr 2 \ r 2 

Here 



2 , 4m 2 c-c : 5A 2 d A 2 t d 2p' . 



^ = / (' lc + ^ + W^ + ^¥ ( c + 7 )1 ' (3 ' 9) 



and the tortoise coordinate r* is defined as 

dr A 
dr* ' c 2 p 2 r 



c (! " TOT ) = f ■ ( 3 - 10 ) 



Notice that for constant p, c one recovers the equations in [27, ^g], as one should. Now, it is quite easy to 
present an example of background flow for which the instability is triggered: take for instance a flow for 
which p is almost constant at infinity (almost means that it asymptotes to a constant value more rapidly 
than the sound velocity). Assume also that, near infinity, c = c\ + £7. Then, we get that near infinity the 
effective potential behaves as 

(3.11) 

r 

For this to have a positive derivative, one requires C2 < [c\ must be positive, as it is the asymptotic 
value of the sound velocity). We thus have one example of flow for which the instability is active. There 
are many others, of course, and there are also instances for which the system is stable. But the important 
point here, is that we can build this effective geometry in the lab, and thus observe this instability. 



- 13 - 



4. The endpoint of the instability 



We have shown that many rotating black objects are unstable against the instability displayed here. We 
have also shown that, for a fixed mass term fj,, the instability disappears for a rotation parameter a below 
a certain value, for which the superradiant factor Re [a;] — J^'i is zero. Now, as the unstable mode 

grows, rotational energy is being extracted from the black brane, and thus its rotation decreases (the 
same phenomenon happens for the black hole bomb §). This means that eventually the instability stops 
growing, when the rotation is in the critical point, for which the superradiant factor is zero. It looks very 
tempting and reasonable to assume that the system will settle down to a low-rotation black brane plus some 
radiation around it. Eventually, all this radiation will escape to infinity. The total angular momentum 
of the system is conserved, so this radiation carries the extra angular momentum. This is also clear from 
the fact that only co-rotating modes are unstable. This is one possibility. But bearing in mind that these 
are bound modes, it seems plausible that one other thing could happen, before the radiation escapes to 
infinity: since we are dealing with higher dimensional objects, the end product could be an object with 
two or more components for the angular momentum, that is to say, there could be a transfer of angular 
momentum from one dimension to the other. 

5. Conclusions 

We have shown that there is a large class of extended rotating black objects which are unstable against 
massless field perturbations. The instability is caused by superradiant amplification of waves trapped 
inside the potential well, or from a different view-point, it is caused by superradiant amplification of wave 
packets bound by the gravitational force. This mechanism is very similar to other superradiant mechanisms 
in the presence of black holes ||, ||, |§, or other rotating objects such as stars p9|1 . The important point 
here is that the extra dimensions work as an effective mass for any massless field. Thus, the evolution of 
a massless scalar field in, for example, the black brane ([l]^) is equivalent to the evolution of a massive 
field in the background of a Myers-Perry rotating black hole. In four dimensions, the black hole displays 
superradiance and has bound states due to the mass term, thus the geometry is unstable. In higher 
dimensions, the spacetimes dealt with here have no stable circular geodesies (the mass term does not give 
rise to a potential well in the wave description) and therefore these spacetimes are stable. This simple 
reasoning implies that general four dimensional extended rotating black objects (more general than the 
one described by (|1,1|)) will be unstable. It also seems to imply that general higher dimensional rotating 
black objects are stable. Moreover, following Zel'dovich 0, it is known that not only the Kerr geometry, 
but any rotating absorbing body for that matter displays superradiance. Thus, this instability appears 
also in analogue black hole models, as we have shown. 

Here, we have derived the instability timescale for a scalar field, and not for geometry (metric) pertur- 
bations, since Teukolsky's formalism for higher dimensional rotating objects is still not available. Still, the 
argument presented above makes it clear that the instability should be present for metric perturbations 
as well. Indeed, the presence of the potential well, in four dimensions, is due only to the mass term, and 
not to the geometry properties themselves, and since metric modes also scatter superradiantly |^] then 
everything we discussed here translates immediately to metric perturbations. We also expect that the 
instability will be stronger for metric modes, because of the following simple reasoning. Superradiance is 
responsible for this instability, and thus the larger the superradiant effects, the stronger the instability. 
Now, we know that in the four dimensional Kerr geometry scalar fields have a maximum superradiant 
amplification factor of about 2%, whereas gravitational modes have maximum superradiant amplification 
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factor of about 138% |19| . Thus, we expect the instability timescale to be almost two orders of magnitude 
smaller for gravitational modes, which means that the instability is two orders of magnitude more effective 
for metric modes. For general black branes, the Gregory-Laflamme instability seems to be stronger than 
the one displayed here, but it is known that certain extremal solutions should not exhibit the Gregory- 
Laflamme instability |$(|, whereas the instability dealt with here should go all the way to extremality. So 
eventually it takes over the Gregory-Laflamme instability. Moreover, recent studies [31| seem to indicate 
that black strings in a Randall-Sundrum inspired 2-brane model do not exhibit the Gregory-Laflamme 
instability, but they should be unstable against this mechanism. 

It seems plausible to assume that the instability will keep growing until the energy and angular momen- 
tum content of the field approaches that of the black brane, when back-reaction effects become important. 
The rotating brane will then begin to spin down, and gravitational and scalar radiation goes off to infinity 
carrying energy and angular momentum. The system will probably be asymptotic to a static, or very slowly 
rotating, final state consisting of a non-rotating black p-brane and some outgoing radiation at infinity. The 
end product could also be an object with two or more components for the angular momentum, that is to 
say, there could be a transfer of angular momentum from one dimension to the other. 
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A. Factorized potential analysis 

We now present a full description of our search for bound states in general (4 + n)-dimensions. We conclude 
that no unstable bound states exist for n > 5. The massive scalar field equation in 4+n dimensional 
spacetime is given by 

Ar~ n 4- (&r n ^-\ +V R = 0, (A.l) 



dr \ dr r 
where Vq is the effective potential, given by 

Vo = {uo(r 2 + a 2 ) - am} 2 - A L 2 r 2 + A lm - 2muoa + uo 2 a 2 + 3U+n-l)a 2 \ . (A.2) 

For simplicity, let us assume Ai m to be A[ m = 1(1 + 1). Note that this assumption is only valid in the limit 
of a 2 (uj 2 — /U 2 ) — > 0. Introducing a new variable and a normalized function <£, defined as 

2 2 

dn = r dr, R = {(r 2 + a 2 )r n }-h$ l (A.3) 

we can reduce equation ( A.2 ) to 

d 2 § 

— + V$ = 0, (A.4) 
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where 



Since V is, due to the simplified Ai m , a quadratic function of u, it can be factorized as 

V = a(u - V+)(u - V-) , (A.6) 
where V + > V_ and a is a positive function, given by 

( r 2 + fl 2)2 _ Aa 2 



ft 



(r 2 + a 2 ) 2 

( r 2 + a 2)2_ a 2( r 2 + a 2_ r l-n) 

(r 2 + a 2 ) 2 
(r 2 + a 2 )r 2 + aV~ n 



(r 2 + a 2 ) 2 



>0. (A.7) 



Note that a also satisfies 



lim a = lim a = 1 . (A.£ 

r— >r_|_ r— >oo 



The factorized potentials V± have the following properties: 

TTICL 

lim V± = -x x = mQ , lim V± = ±(Jt, (A. 9) 

where f2 is the rotation frequency of the black hole at the horizon. Note that V+ = V_ is satisfied at the 
horizon. 

The factorized potentials are quite useful to see how a solution behaves when the frequency u is given. 
First, let us consider the four dimensional case, i.e. the n=0 case. Taking the case of I = m = 1, j = 0, 
fj, = 0.7, and a = 0.4999 as an example, we plot the typical behavior of the factorized potentials for the 
four dimensional case in Figure ||. Since a is positive definite, the field $ has propagative (evanescent) 
character when oj > V+ or u < V- (VI < u < V+) because ^(f^/drl < (§~ l d 2 $/drl > 0) there. In 
other words, the region where u > V+ or u < V_ (V_ < cj < V+) is a classically allowed (forbidden) one. 
By a factorized potential diagram, we can then see where is the propagative or evanescent zone once the 
frequency u is given. 

Next, let us consider bound states for our problem. To have a bound state, we need a propagative 
zone between evanescent zones. That is to say, unless V + (V_) does have a local minimum (maximum), 
there are no bound state in the potential. We can therefore expect from the behavior of dV±/dr whether 
bound states are possible or not. In Figure ||, we display the dV^/dr = curve on the r _1 - fi plan for 
the case of I = m = 1, j = 0, and a = 0.4999. From the figure, we can confirm that V+ as a function of r 
have two extremes for < \i < 1.321. From Figure ^, we can indeed see that V+, in which \i = 0.7, has 
a local minimum near r _1 = 0.15 and that there is a potential bound state if we take oj ~ 0.6878, which 
is a correct eigenfrequency of the bound states. Another important thing is that u ~ 0.6878 is below m$7, 
which can be easily seen from the factorized potential diagrams. This means that this bound state is an 
unstable mode because it satisfies the condition of the superradiant instability. 
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Figure 4: The factorized potentials V± as a function of r 1 . The potentials are characterized by n = 0, I — m = 1, 
j = 0, fi = 0.7, and a = 0.4999. 




Figure 5: The curve of dV+/dr(r 1 ,fi) — for the parameter n = 0, I = m = 1, j = 0, and a = 0.4999. 



Now, let us move on to higher dimensions, taking n = 1 to be specific. Consider the case of I = m = 1 
and j = as an example. We are interested in finding out whether or not there can be bound states in 
the five dimensional case. In order to see the distribution of extremes of V+, which give us a criterion for 
nonexistence of the bound states, a curve of dV+/dr(r, fj,) = on the r~ l - fj, plan is sketched in Figure^, 
where the black hole characterized by a = 0.9 is taken. In the figure, it is found that if the mass \x satisfies 
1.658 < fi < 1.784, dV+/dr as a function of r has two zeros; one corresponds to a local maximum of V+ 
and the other a local minimum of V+. It should be emphasized that contrary to the four dimensional case, 
there is a critical mass below which no potential bound state exists. Therefore, there is the possibility 
that a bound state exists in the five dimensional case. A distribution of V+ for the case of \i = 1.75 is 
shown in Figure ^. We can conclude from the figure that all the frequencies for potential bound states are 
always above mO. Therefore, those potential bound states must not be unstable against the superradiance 
instability even if they exist. For other sets of the parameters, we have found similar behavior of the 
factorized potential. This means that in the five dimensional case, there might be bound states but all the 
potential bound stats are stable. This analysis is consistent with our numerical studies, in which we have 
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Figure 6: The curve of dV+/dr(r~ x , fi) = for the parameter n = 1, I = m = 1, j = 0, and a = 0.9. 
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Figure 7: The factorized potentials V± as a function of r _1 . The potentials are characterized by n= 1, I = m = 1, 
j = 0, fi = 1.75, and a = 0.9 (left panel). The right panel is the same as the left one, but, a close-up of the potential 
V+ near extremes. 

not found any unstable modes. 

For general n > 5 the trend is the same, there are no unstable bound states. 
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